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Berry and Wang [Phys. Rev. A 83, 042317 (2011)] show numerically that a discrete-time quan- 
tum random walk of two noninteracting particles is able to distinguish some non- isomorphic strongly 
regular graphs from the same family. Here we analytically demonstrate how it is possible for these 
walks to distinguish such graphs, while continuous-time quantum walks of two noninteracting parti- 
cles cannot. We show analytically and numerically that even single-particle discrete-time quantum 
random walks can distinguish some strongly regular graphs, though not as many as two-particle 
noninteracting discrete-time walks. Additionally, we demonstrate how, given the same quantum 
random walk, subtle differences in the graph certificate construction algorithm can nontrivially im- 
pact the walk's distinguishing power. We also show that no continuous-time walk of a fixed number 
of particles can distinguish all strongly regular graphs when used in conjunction with any of the 
graph certificates we consider. We extend this constraint to discrete-time walks of fixed numbers of 
noninteracting particles for one kind of graph certificate; it remains an open question as to whether 
or not this constraint applies to the other graph certificates we consider. 



I. INTRODUCTION 

The graph isomorphism problem is considered a good 
candidate for speedup by a quantum algorithm. The 
graph isomorphism (GI) problem is, given two graphs, 
to determine if they are isomorphic, that is, if one can 
be mapped to the other via a relabeling of vertices. To 
date, the best known classical algorithm has a runtime of 
0(c^ l ° sN ), where c is a constant and N is the number 
of vertices in the graph [1 . GI is thought to be simi- 
lar to integer factoring [2] in that both could reside in 
the complexity class NP-Intermediate. Given that there 
exists a polynomial-time quantum algorithm for integer 
factoring, it is thought that there may exist a quantum 
speedup for GI as well. 

One class of algorithms that has been explored for GI 
is that of quantum random walks. Quantum random 
walks (QRWs) are of particular algorithmic interest as 
they are known to offer certain computational advantages 
over classical random walks [3H5]. There exist problems 
for which QRW algorithms are known to have faster run- 
times than their fastest known classical analogues [6HI2] . 

Much work has been devoted to developing a quan- 
tum algorithm for the graph isomorphism problem. Sev- 
eral works have examined ways to use continuous-time 
quantum random walks (CTQRWs) to solve GI [T3lU5j . 
Others have looked to utilize discrete-time quantum ran- 
dom walks (DTQRWs) for GI algorithms [TBI Q7]. While 
these efforts have helped illuminate what kinds of QRW 
algorithms might be fruitful for GI, an efficient QRW al- 
gorithm for GI has yet to be developed. Moreover, it is 
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unknown if DTQRWs or CTQRWs are better candidates 
for an efficient GI algorithm, or if both classes of walks 
will ultimately be of equal utility. 

This question is explored in a recent paper by Berry 
and Wang [TB]. Their results indicate that discrete-time 
quantum random walks may have greater ability to dis- 
tinguish non-isomorphic strongly regular graphs (a class 
of graphs which are particularly difficult to distinguish) 
than continuous-time quantum random walks. This is 
based on (1) the ability of a specific discrete-time quan- 
tum random walk of two noninteracting particles to dis- 
tinguish strongly regular graphs (SRGs) that are not dis- 
tinguished by an algorithm based on a continuous-time 
quantum random walk of two noninteracting particles 
|14j . and (2) the fact that the dimension of the state 
space of discrete-time walks on graphs is larger than that 
of continuous-time walks. 

In this paper, we explain why the two-particle non- 
interacting walks of Berry and Wang are able to dis- 
tinguish non-isomorphic SRGs. To do this, we show 
analytically how even single-particle DTQRWs are able 
to distinguish non-isomorphic SRGs. This demonstrates 
how two-particle noninteracting DTQRWs can also dis- 
tinguish SRGs. We then verify numerically that this dis- 
tinguishing power of single-particle DTQRWs exists. 

Additionally, we explore the differences between the 
comparison algorithm used by Berry and Wang and the 
one used by Gamble et al. [14] . Both comparison algo- 
rithms rely on the construction of sorted lists of data, or 
"graph certificates" , which are used to compare graphs. 
Though the respective certificates are constructed in sim- 
ilar manners, we find that, when applied to single-particle 
discrete-time quantum random walks, there is a non- 
trivial difference between the distinguishing power of the 
graph certificates of [14] and [18]. Thus we demonstrate 
the importance of using the same kind of graph cer- 
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tificates when attempting to compare the distinguishing 
power of DTQRWs and CTQRWs. 

Lastly, we extend a result of Rudinger et al. [H] , who 
showed that there does not exist a CTQRW with a fixed 
number of noninteracting particles that can distinguish 
all SRGs when used in conjunction with the comparison 
algorithm of Gamble et al. Here we show that this ar- 
gument holds for DTQRWs as well; no DTQRW with a 
fixed number of noninteracting particles can distinguish 
all SRGs when the comparison algorithm of Gamble et 
al. is used. We also examine the comparison algorithm 
of Berry and Wang, as well as an algorithm with similar- 
ities to both the algorithms of Berry and Wang as well as 
Gamble et al. We find that no CTQRW of a fixed num- 
ber of noninteracting particles can distinguish all SRGs 
when used with these comparison algorithms. It remains 
an open question as to whether or not this is true for 
DTQRWs of a fixed number of noninteracting particles. 

This paper is organized as follows: Section [n] pro- 
vides the requisite background regarding graph isomor- 
phism, strongly regular graphs, and quantum random 
walks. Section |III| describes in detail the graph compar- 
ison procedures of Gamble et al., and Berry and Wang, 
as well as a procedure which is "in between" the two. In 
Section [TV] we analytically show how the single-particle 
DTQRW has the potential to distinguish non-isomorphic 
graphs when used with any of the three comparison algo- 
rithms we consider. Section |V| provides numerical results 
for single-particle DTQRWs; we find that all three com- 
parison algorithms numerically distinguish some, but not 
all, of the tested SRG pairs, and one comparison proce- 
dure is not as strong as the other two. In Section VI 



we show that a DTQRW with a fixed number of non- 
interacting particles cannot distinguish all SRGs when 
used with the comparison procedure of Gamble et al.; we 
show the same result holds when the comparison algo- 
rithm of Berry and Wang is used with CTQRWs with a 
fixed number of noninteracting particles. We are unable 
to extend this result for Berry and Wang's comparison 
algorithm when used with a DTQRW of a fixed number 
of noninteracting particles, indicating a potential differ- 
ence between DTQRWs and CTQRWs. We discuss our 
results in Section |VII| Lastly, the Appendix provides the 
explicit formula for a quantity necessary to demonstrate 
the distinguishing power of the comparison algorithm of 
Berry and Wang. 



unordered pairs of vertices. Two vertices are connected 
by an edge, and are said to be adjacent or neighboring, if 
and only if that vertex pair appears in E. A convenient 
representation of a graph is its adjacency matrix A: 



1 if vertices i and j are adjacent 
An = < (1) 
if vertices i and j are not adjacent. 



A graph of N vertices has an N x N adjacency matrix. 
As all graphs we consider are simple and undirected, here 
A is always symmetric, with zeros on the diagonal. 

Two graphs are isomorphic if and only if there exists 
a relabeling of vertices which transforms one graph into 
the other. More formally, two graphs represented by ad- 
jacency matrices A and B respectively are said to be iso- 
morphic if and only if there exists a permutation matrix 
P such that B = P^AP. 



B. Strongly regular graphs 

The graphs considered by Berry and Wang, as well 
as Gamble et al. and Rudinger et al. are called strongly 
regular graphs (SRGs). Strongly regular graphs are good 
test cases for any candidate algorithm for graph isomor- 
phism, as there exist many relatively small SRGs which 
arc difficult to distinguish with the best known classical 
algorithms. 

SRGs are described by four parameters, denoted 
(N,d,X,fj,). (Sometimes k is used instead of d.) N is 
the number of vertices in the graph, and every vertex 
is adjacent to d other vertices (the graph is d-regular). 
Every two adjacent vertices share A common neighbors, 
while every pair of non-adjacent vertices share /i com- 
mon neighbors. SRGs with the same parameters are said 
to be in the same SRG family. Infinitely many SRG 
families have more than one non-isomorphic graph [19] . 
Since family parameters are easily computed, it is non- 
isomorphic graphs in the same family that are difficult 
to distinguish. SRG families that contain multiple non- 
isomorphic graphs are studied in P31 [T5J HH] • 

It follows from the definition of SRGs that SRG adja- 
cency matrices satisfy 20J : 



(d - fi)l + /iJ + (A - fx) A, 



(2) 



II. BACKGROUND 

A. Basic graph definitions 

Here we review the graph-theoretic background neces- 
sary to discuss quantum random walks on graphs. As in 
[IH [15l [18] , we consider only simple, undirected graphs. 
A graph G — (V, E) is a set of vertices V and edges be- 
tween those vertices E. The set of vertices V is a set 
of labels, usually integers, and the edge set E is a set of 



where 1 is the identity matrix and J is the matrix of all 
ones. Given that J 2 = NJ, 3 A = AJ = dA, and 1 acts 
trivially on 1, J, and A, we see that {1,J,A} forms a 
commutative three-dimensional algebra. Thus, for any 
non-negative integer n: 



A™ = a n t + (3 n 3 + 7 „A, 



(3) 



where a n , f3 n , and 7„ all depend only on the SRG family 
parameters and n. 
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C. Discrete-time and continuous-time quantum 
random walks of noninteracting particles 

Next we discuss how to form quantum random walks 
on graphs. For the noninteracting continuous-time model 
used by [HI US] , a Hamiltonian is defined for a graph of 
N vertices with adjacency matrix A: 

N 

H=-^A <J -cJc J , (4) 

where e] and Cj are the creation and annihilation oper- 
ators, respectively, for a boson or (spinless) fermion on 
site i. For bosons, they satisfy the commutation rela- 
tions [cj, c^} = 6ij and [cj, Cj] — [cj, cj] = 0. For fermions, 
they satisfy the anti-commutation relations {ci,c|} = Sij 

and {ci,Cj} = {c|,ct} = 0. This Hamiltonian is a tight- 
binding model, where each site corresponds to a vertex 
in the graph; a particle may move from one vertex to an- 
other if the vertices are adjacent. This Hamiltonian may 
be used with any number of particles (either bosons or 
fermions). A matrix representation of this Hamiltonian 
for a fixed number of particles is given in |15j . The evo- 
lution operator is then defined in the standard fashion: 

U c (t) = c-' itH , (5) 

where we have h = 1 for convenience. The subscript C 
denotes that the evolution operator is for a continuous- 
time walk; discrete-time evolution operators will be de- 
noted with a subscript D. 

For discrete-time walks, each particle has an auxiliary 
coin state, which "points" to where the particle will next 
move. Thus the single-particle walk has basis states of 
the form \ij), denoting the particle is on site i, and its 
coin is on site j, which by definition must be adjacent to 
i. 

The noninteracting walk considered by 18J is based 
on a single-particle Grover-coined walk. We follow [TS] 
and express the evolution operator in expanded form, 
allowing the coin index to run over all vertices in the 
graph, but still requiring the particle be adjacent to its 
coin. Thus, for an TV- vertex graph with adjacency matrix 
A, the single-particle discrete-time evolution operator is 
given by 

U D = SA(1 (g) C)A. (6) 

A\ij) = Aij ensuring that the coin is adjacent to 
the particle. S is the swap operator; S \ij) — C is 
the coin operator. Here it is the Grover coin, which for 
a d- regular graph is C = —1+ |j. This expanded form 
of Ufl has the same behavior as in [181 ; we have only 
introduced N 2 - Nd rows of 0, so its dimension is now 
N 2 x N 2 [T5]. 

To advance the walk a discrete number of time steps 
t, XJd is applied t times: Uo(t) = (Uo)'. Additionally, 
the evolution operator for p noninteracting particles is 
just p tensor copies of U^: Ud, p — (U_d)® p . 



III. COMPARISON ALGORITHMS 

Now that we have defined evolution operators for 
continuous-time and discrete-time walks, we explore how 
these operators can be used to compare graphs and test 
for isomorphism. All methods of comparing graphs that 
we examine are based on, given an evolution operator, 
generating a list of numbers sorted by size for each of the 
graphs, and comparing the respective lists. We follow 
[IB] and refer to these sorted lists as graph certificates. 

There are multiple graph certificates that can be con- 
structed from a given evolution operator. In this sec- 
tion, we examine three classes of graph certificates. First, 
we give the definition of the certificate which is used by 
[HI [15] for CTQRWs. We then describe a certificate 
which is "between" the certificate of [T31 [IS] and the cer- 
tificate of [IS]. Lastly, we define the certificate used by 
[T%] for DTQRWs. All three certificates can be applied 
to either continuous-time or discrete-time QRWs. 

As the certificates are lists, we denote the three dif- 
ferent certificate classes as Lq, L\, and L2, respec- 
tively. To indicate that a specific certificate corresponds 
to a discrete-time walk, a superscript D is included; a 
superscript C denotes the certificate corresponds to a 
continuous-time walk. Additionally, a second superscript 
is used to indicate that number of particles in the walk 
which generates the certificate. For example, L®'^ would 
refer to the Li graph certificate for a discrete-time walk 
of a single particle. 

We call the graph certificate of [T3] and [TS] Lq. It is 
defined as follows: 

L Q (A,t) = sort ({|Z7(i) mn | : Vm, n e dimU}) . (7) 

Thus, if the possible values the evolution operator ele- 
ments can take on (up to different phases) are different 
for the different graphs, then the graphs will be distin- 
guished. Even if the possible values each element can take 
on are the same for the different graphs, the graphs will 
be distinguished if those values have different multiplici- 
ties in the two evolution operators. While this algorithm 
has been used for CTQRWs, it can in principle be used 
for DTQRWs as well. Additionally, we see from its def- 
inition that Lq is naturally defined for a walk with any 
number of particles. 

O 1 

It was proven in [13] that L ' certificates could not 
be used to distinguish non-isomorphic SRGs from the 
same family, while [14] extended this proof for CTQRWs 
of two noninteracting particles. However, |15j showed 
that Lq for CTQRWs with three or more noninteracting 
particles could distinguish many (but not all) pairs of 
non isomorphic SRGs from the same family. 

The next graph certificate is designed for DTQRWs, 
but can be used for CTQRWs. We denote it Lx, which 
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for the single-particle DTQRW is defined as 



/v 



Lf'\A,t) = sort J2 I fe'l U AiW \kl) 
Vi,k,l€{l...N}\ 



(8) 



Each element of if' 1 represents the total probability of 
a particle being on vertex i after t steps, given an initial 
state \kl) . A natural extension of this method for a two- 
particle CTQRW is 



Therefore, their graph certificates (which we denote 

L2 ) have at most nonzero elements, as a d-regular 
graph with TV vertices contains edges. The certificate 
is defined as 

~D.2 



ixhiih\ U A2 (t)^(|Mifc) + \lkkl))f 
" ' (13) 




'=1 jlj2 = l 

Vii,i 2 ,G {l...JV},V(*,/) G E 



Lf' 2 (A,t) = sort 




Vi,fc € {l.../V},V7 e {l...fc} 



(9) 



where E is the edge set of A. The physical interpretation 

~d,2 n „ 

of the certificate L 2 is identical to that of L 2 : except 

that each initial state is delocalized across an edge in the 

graph [18]. 



where \ij) B denotes the bosonic (symmetrized) state in 
which there is a boson on site k and a boson on site I. 
(This method may be used just as well with fermions, 
with the basis states appropriately anti-symmetrized.) 
The extension of L\ for more than two continuous-time 
particles is given in Section [VlB| 

The final graph certificate we examine is a variant 
of this method, and is the one utilized by |18| . This 
certificate, denoted L2, is defined for the single-particle 
DTQRW as follows: 

L°' 1 (A,T) = 

sort ( I E E 1 (y 1 u ai (*) 1 kl ) 1 2 : v *> k > 1 e {i • ■ ■ n} I y 

(10) 

This method sums up the probabilities of a particle being 
at a particular site at different times, given the same 
initial state. For the two-boson CTQRW, this certificate 
has the form: 



N 



L^ 2 (A,T) = sort! 



t=i j=i 



Vi,k e {l...JV},Vi €{!...*} 



(11) 



IV. ANALYTIC DEMONSTRATION OF THE 
DISTINGUISHING POWER OF DTQRWS ON 
SRGS 

Berry and Wang demonstrated numerically that the 
noninteracting two-particle DTQRW with the L 2 certifi- 
cate method could distinguish many SRGs. In this sec- 
tion, we show analytically why this is possible. To do 
so, we first analytically explore the distinguishing power 
available to single-particle walks using only the L and 
Li graph certificates. 

In general, if a particular kind of certificate will always 
fail to distinguish two non-isomorphic SRGs, it is because 
all elements of a certificate, as well as their multiplicities, 
are functions of SRG family parameters. This is how we 
show, or fail to show, the limitations of each certificate 
considered in this section. 



Single-particle DTQRW with L graph 
certificate 



Here we show analytically that even single-particle 
DTQRWs have the potential to distinguish SRGs, us- 
ing only the L graph certificate. To begin our analysis, 
we compute an arbitrary element of XJd.i for t = 1: 



As with L\i L2 can be extended to any number of 
discrete-time or continuous-time walkers. 

For a two-particle DTQRW, the size of L 2 ' can be 
quite large, with as many as N 4 d 2 non-zero elements. In 
[15] the initial states are limited to what Berry and Wang 
call "bosonic edge states" . Given two adjacent vertices k 
and I, a bosonic edge state, denoted by Berry and Wang 
as \(3 + ), is defined as 



1 



V2 



(\kllk) + \lkkl)). 



(12) 



Ud,i \kl) = A i:j A k id jk 



(14) 



if Aij = 
because 



The particle on site k with 
move to site i with its coin 

Aki = Sjk = 1- 
both the bra and 
state must always have its 
adjacent to the location of 
of this transition is equal 
otherwise. We see that the 



its coin pointing to site I can 
pointing to site j if and only 
(Aij = A^i = 1 is required 
ket must be legal states; a 
coin point to a site that is 
the particle.) The amplitude 
to I if i = I, and -1 + | 
particle's movement for this 
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single time step is highly restricted, and that all possible 
evolution operator element values are strictly functions 
of the family parameters. 

Now that we have computed the possible non-zero val- 
ues for U \ kl), we compute the multiplicities of these 
values. If the multiplicities of one value are different for 
two SRGs in the same SRG family, then the Lq' 1 certifi- 
cates will distinguish the graphs. 

To compute the multiplicities, we follow []3J[T5]. We 
denote by M(x) the multiplicity of the value x in the evo- 
lution operator, and we recall that the family parameters 
for an SRG are denoted (N, d, A, /j,). We find: 



When all four vertices are connected to each other, 



M 



AI 



N 



AijA k i6jk(l 

ijkl 



Su), 



N 

AijAkiSjkSu. 

ijkl 



(15) 



(16) 



Each of these summands is a product of four terms, 
each of which corresponds to an identity or adjacency 
relationship that appears in \5d,i \kl). 1 — Su ap- 
pears in the first summand, because Su is "turned off" 
(equals 0) when (ij | U^^ \kl) = §; Su appears in the 
second summand because Su is "turned on" (equals 1) 
when Ud,i \kl) = -1 +_§. 



The sums given in Eqs. ( 15 ) and ( 16 1 can be computed 



straightforwardly because all the index contractions here 
are reducible to matrix multiplication and traces. Addi- 
tionally, we use the SRG identity of Eq. ((2), and find 



M 



— Nfi(N — d—1) + NXd, 



M -1 



= Nd. 



(17) 



(18) 



We see that these multiplicities are functions of the fam- 
ily parameters, so the walk has no distinguishing power 
when t = 1. 

Next, we consider later times. U^, l \kl) for t = 2 
and t — 3 are straightforward to calculate; both the val- 
ues and degeneracies for these cases are reducible to sums 
over products of adjacency matrices and traces of adja- 
cency matrices, and therefore can be written as func- 
tions of the SRG family parameters. Therefore, \5 2 D 1 
and U|, l cannot distinguish non-isomorphic SRGs of the 
same family when used to generate Lq certificates. 

At t = 4, all six possible adjacency relations appear 
(Aij, Am, A jk , Aik, Aji and An): 

(ij\U 4 \kl) = A t] A kl {4d- 2 (A u - A jk ) + 2cT 1 x (19) 
(2S ik - A u S lk - A jk 5 ik + Sji - AiiSji) + S ik Sji~ 
8d~ 3 ((d - n)(Sik + Sji) + (A - fi)(A ik + A jt )+ 
2fi) + 16d- 4 (^- fc (d - n)(\ - im)+ 
A jk (d + (A - M ) 2 - fi) + (d + A - p)n)). 



\kl) 



-16A 



+ 



w{d+(\- fJ .y- t j.+(d+\- tJ .) t i) 



As 



df 1 d 4 
no other configuration of vertices yields this value for 

\ U|j 1 \kl), this value appears in the operator \J A D 1 M 

times, where M satisfies 



N 

M = AijAkiA^AtiAjkAj 

ijkl 



(20) 



Computing the sum in Eq. (20) requires contracting 
over four indices, each of which occurs three times. Such 
a sum cannot be reduced to be in terms of matrix mul- 
tiplication and traces. Ref. |15) showed, in the context 
of continuous-time noninteracting walks of three parti- 
cles, that these kinds of sums are functions of the num- 
ber of shared neighbors belonging to triples of vertices in 
the graph in question. This number is, in general, dif- 
ferent for different triples of vertices in the same SRG, 
and therefore not a function of the SRG family param- 
eters. Therefore, XJj-, x has the potential to distinguish 
non-isomorphic SRGs when used with the Lq method. In 
Section |Vj we numerically show that this method does 
indeed distinguish some (but not all) non-isomorphic 
SRGs. 

In contrast, in Refs. (T3l [14] it is proved that single- 
particle and noninteracting two-particle CTQRWs with 
the Lq graph certificate method cannot distinguish any 
non-isomorphic SRG pair from the same family. Sev- 
eral key differences between DTQRWs and CTQRWs be- 
come apparent. First, an element of the single-particle 
continuous-time walk evolution operator is indexed by up 
to 2 vertices (as the particle is on 1 vertex in an initial 
state, 1 vertex in the final state, and the sets of final and 
initial vertices need not overlap) . However, an element of 
the corresponding single-particle discrete-time walk evo- 
lution operator is indexed by up to 4 vertices, because 
each particle corresponds to a vertex, as does each par- 
ticle's coin. Additionally, in the continuous-time walk, 
the value of an evolution operator element does not de- 
pend on whether or not two vertices which are both in 
the final or initial state are adjacent. Therefore, such 
adjacency relations are not considered when performing 
the appropriate multiplicity sum. However, because a 
particle must always be adjacent to its coin, the adja- 
cency relation between a particle and its coin is always 
included in the discrete-time sum used to compute cle- 
ment multiplicity. The presence of these additional ad- 
jacency relations gives the single-particle DTQRW with 
the Lq comparison protocol the potential to distinguish 
non-isomorphic SRGs. 



B. Single-particle DTQRW with Li graph 
certificate 



Here we show that the Li method, when used with the 
single-particle DTQRW, also has the potential to distin- 



6 



guish non-isomorphic SRGs from the same family. To inspection of Eqs. (18]) and (10), we see that for a final 



start, we compute an element of Y^=i I 1 Ud,i \kl) 



for t = 1. 



N 

E 



<y| U Dl i |fc/> | 2 = A lk A kl (8 a (l - f ) + £) (21) 



Possible non-zero values for SjLi I Ud,i \kl) | 2 are 



r/- 



and 1 



which are both functions of the family 



parameters. We may compute their multiplicities, and 
find that they are also functions of the family parameters. 
Thus, L( ) ' 1 (< = 1) graph certificates cannot distinguish 
non-isomorphic SRGs from the same family. 

Similarly, it may be shown that at t = 2 for the single- 
particle DTQRW, all elements of L\ are functions of the 
family parameters, as well as their multiplicities. There- 
fore, if' 1 (i = 2) certificates also cannot distinguish non- 
isomorphic SRGs from the same family. 

However, at t = 3, the values of the elements of Li for 
the single-particle DTQRW at t = 3 are not functions of 
the family parameters. The Appendix provides the value 
of J^7=i I Ui? i i(3) \kl) | 2 , which has the form: 



N 

E 



N 

U A i(3) \kl) | 2 = g(i, k, l)+h(i, k, I) AuAflAjk, 

i=i 

where g(i,k,l) and h(i,k,l) are functions of the fam- 
ily parameters. However, X!j=i AijAjiAjk is not a 
function strictly of family parameters, for the same 
reason that Eq 
need to examine t 



(20) 



not. Thus we do not even 
ic multiplicities of different values of 
Yl'j=i I Ud,i(3) \kl) | 2 , as this sum takes on values 
that are not functions of the SRG family parameters. 
Therefore, the single-particle L\ method at t = 3 can 
potentially distinguish non-isomorphic SRGs from the 
same family. In Section [Vj we numerically show that 
this method can distinguish some non-isomorphic graph 
pairs from the same SRG family. 

However, this is not true for two-particle noninteract- 
ing CTQRWs. Using the methods of [H EDS, it can be 
shown that, when the appropriately symmetrized or anti- 
symmetrized states are used with the continuous-time 
evolution operator, all possible values and corresponding 
multiplicities of X)j=i I (*J I Uc,2(i) \kl) | 2 are functions of 
the family parameters and the time t. Thus, we see 
that the single-particle DTQRW L\ method is strictly 
stronger than the two-particle noninteracting CTQRW 
Li method for distinguishing SRGs. 



C. Single-particle DTQRW with L 2 graph 
certificate 

Now we examine L^ > ' 1 (T) analytically for varying val- 
ues of T. Recall that elements of L^^lT) are sums of 
elements from L^' 1 (i) with t running from 1 to T. By 



time of T = 1, L 2 is the same as L\(t = 1), and will 
have no distinguishing power. Similarly, for T = 2, each 
element of L 2 will be a sum of a term from L\(t = 1) 
and L\(t = 2). As the value and multiplicity of each 
of those terms is strictly a function of SRG family pa- 
rameters, the corresponding values and multiplicities of 
each element of L2(T = 2) will be a function of SRG 
family parameters. Thus L^ ),1 (T = 2) certificates cannot 
distinguish non-isomorphic SRGs from the same family. 

However, we know that L® :1 (t = 3) certificates have 
the potential to distinguish SRGs. Thus, L^ >,1 (T = 3) 
certificates also have the potential to distinguish SRGs 
from the same family, as the elements of a L^ > ' 1 (T = 3) 
certificate are sums of terms which include elements of a 
L^' 1 ^ = 3) certificate, which we have shown to not be 
functions of family parameters. 

Two-particle noninteracting CTQRWs with the L 2 
method may be similarly analyzed. Because the L\ 
method fails to distinguish SRGs with noninteracting 
two-particle CTQRWs for all values of t, these walks 
with the L 2 method will also be unable to distinguish 
SRGs. Thus we analytically see that for all three compar- 
ison methods contemplated, the single-particle DTQRW 
has the potential to distinguish non-isomorphic SRGs 
from the same family, while two-particle noninteracting 
CTQRWs cannot. 



V. NUMERICAL RESULTS 

We now test the three comparison methods with the 
single-particle DTQRW on three families of SRGs. For 
the DTQRWs, all three algorithms are tested at the min- 
imum time they can potentially distinguish SRGs (t = 4 



for Lq' 1 , t — 3 for L^' x , and T = 3 for 



D.l 



-DM 



We ad- 



ditionally test Lq' 1 at t = 2N, L®' 1 at t = 4, and L®' 1 
at T = 4 and T = 2N. (T = 2N is the total time used 
by Berry and Wang for their tests [IB].) Our results are 
given in Table [I] 

Each algorithm distinguishes some, but not all, non- 
isomorphic SRGs from the same family. Also, no single- 
particle DTQRW algorithm that we test is able to distin- 
guish the pair of SRGs in (16, 6, 2, 2), which is in contrast 
to the results of Berry and Wang, who successfully dis- 
tinguished these graphs using the two-particle DTQRW 

~D,2 

with L2 certificates. Additionally, the various algo- 
rithms we test have significant distinguishing power on 
the other two SRG families we examine. We note that 
the distinguishing power for each algorithm saturates 
at the minimum time required to allow for potentially 
distinguishing SRGs. This could be because SRGs are 
distance-regular with diameter 2, so allowing the parti- 
cles to walk for longer does not actually result in captur- 
ing more of the graph structure. 

The distinguishing power seems to saturate with 
L^ >,1 (t = 3), for higher times, either with L^' 1 or 
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Number of Undistinguished Pairs 


SRG Family 
(N ,d,X,n) 


Number 
of Graphs 


Comparisons 


L^(i = 4) 


L%> 1 (t=2N) 


L?' x (t = 3) 


Lf' 1 (t = 4) 


L°>\T = 3) 


L»>\T = 4) 


L°' 1 {T = 2N) 


(16,6,2,2) 


2 


1 


1 


1 


1 


1 


1 


1 


1 


(25,12,5,6) 


15 


105 


11 


11 

















(26,10,3,4) 


10 


45 


3 


3 


1 


1 


1 


1 


1 



TABLE I. Numerical results for graph isomorphism testing using single-particle DTQRWs with varying comparison algorithms 
with various times. The different graph certificates {L®' 1 , Lf' 1 , and are defined in Eqs. Q, B, and (10 1. We note that 

all three algorithms have significant, but not universal, distinguishing power on SRGs. Additionally, this distinguishing power 
saturates at the minimum time at which each algorithm can potentially distinguish SRGs. Lastly, we see that L 1 ^ and L^' 1 
are more powerful than L ' , but L 2 ' is possibly no more powerful than L 1 ' . 



L 2 ' , no further distinguishing power is obtained. For 
(25, 12, 5, 6), L,f ' l (t = 3) distinguishes all possible graph 
pairs. For (26, 10, 3, 4), there is one graph pair that can- 
not be distinguished by L l ' (t = 3), L x ' (t = 4), or L 2 ,x 
with any time we test. Therefore, it is possible that L^' 1 
offers no more distinguishing power than L®' 1 . More 
generally, it is possible that L 2 ' p yields no more distin- 
guishing power than L®' p . Additionally, it is possible 
that there exists a certain time (and perhaps relatively 
small time) beyond which no additional distinguishing 
power is obtained. As the complexity of computing L 1 ,p 
and L 2 ' p increase with the number time steps used, it 
could useful to know if there is no additional information 
to be gained by increasing the number of time steps. 



VI. ASYMPTOTIC LIMITS FOR DTQRWS AND 
CTQRWS 

We have shown how, by all three comparison algo- 
rithms considered, the single-particle DTQRW has more 
distinguishing power on SRGs than noninteracting two- 
particle CTQRWs (which cannot distinguish any non- 
isomorphic pair of SRGs from the same family). How- 
ever, as the distinguishing power of the single-particle 
DTQRW is not universal on SRGs, we can contemplate 
how to increase this distinguishing power. One obvious 
method is to add additional particles to the walk. In- 
deed, with a two-particle noninteracting DTQRW, Berry 
and Wang distinguished SRG pairs that are not distin- 
guished by L2' 1 certificates for T = 3,4 and 2N [15] , 
Thus it is interesting to characterize the asymptotic dis- 
tinguishing power of many-body QRWs. In this section, 
we show that no CTQRW of a fixed number of noninter- 
acting particles can distinguish all SRGs when used with 
either Lq, Li, or L2 certificates. We also show this result 
holds for DTQRWs of any fixed number of noninteract- 
ing particles when used to generate Lq certificates; we 
are unable to extend this result for such walks with L\ 
and Li certificates. 



A. Asymptotic behavior for Lo 

It was demonstrated in [15] that increasing the num- 
ber of noninteracting particles in a CTQRW could sig- 
nificantly increase distinguishing power of the walk when 
used with the Lq method. However, it was also shown 
in [T5] that no CTQRW with a fixed number of nonin- 
teracting particles could distinguish all SRGs. This was 
done by showing that, as graph size increases, the num- 
ber of unique SRGs with the same family parameters is 
super-exponentially large in graph size |19l I21j , while for 
a fixed particle number p, the number of unique graphs 
that the p-particle CTQRW can distinguish with the Lq 
method is polynomial in graph size. 

We show here that this limitation of Lq certificates 
applies to DTQRWs as well. Let us consider the 
p-particle noninteracting DTQRW, and let us denote 
by X p the maximum number of unique values that 
(ii,ji.--ip,jp\U t DiP \ki,li...kp,lp) can take on. X p is 
not a function of the size of the graph; the value of a 
DTQRW evolution operator element for an SRG is de- 
termined only by the configuration of up to 4p vertices 
in that SRG. (For further details, see [15].) 

Additionally, the number of elements in Ud iP is N ip , 
where N is the number of vertices in the graph. Thus, 
for an SRG family with graphs of N vertices, the Lq 
method generates certificates of length N 4p , where each 
element in the certificate can take on a maximum num- 
ber of X p different values. Thus, given an SRG family 
with graph size N, the maximum number of unique Lq 
graph certificates the noninteracting p-particle DTQRW 
can generate is equivalent to the number of ways one can 
put N Ap indistinguishable balls into X p distinguishable 

bins, or r^-!- 1 ) E2- 

As Xp is constant in TV, this quantity is polynomial in 
N. Thus, just like noninteracting CTQRWs with fixed 
particle number, DTQRWs with fixed particle number 
cannot distinguish all SRGs with the Lq graph certifi- 
cate method. Thus, the Lq method for noninteracting 
p-particle walks cannot yield a universal GI algorithm 
for SRGs, whether or not the walks are continuous-time 
or discrete-time. 



B. Asymptotic behavior for L\ and L2 

For discrete-time walks, the L\ graph certificate does 
not necessarily have the same limitations as L . This 
is because the values of the elements of discrete-time L\ 
certificates are not functions of SRG family parameters, 
as demonstrated by Eq. (22). For a fixed number of 



discrete-time walkers, the number of unique values that 

, ,, , , , 1 t jt 



E Jx ...j I I ■ ■ ■ ip3p\ V D,p I Ml 



kplp) 1 2 can take on 



is not a function only of particle number. In principle, 
this quantity may be a function of the graph size, in 
which case the number of unique L\ graph certificates for 
a fixed particle number and SRG family may be super- 
exponential in graph size. Thus, the number of unique 
SRGs in a family with graph size N is not guaranteed 
to be larger than the number of unique discrete-time L\ 
graph certificates that can be generated for walks of p 
nonintcracting particles on graphs of TV vertices. Thus, 
the proof of the limitations of the L method fails to 
translate to the L\ method for DTQRWs. Whether there 
exists such a limitation is therefore still an open question. 

Because DTQRW L 2 graph certificates are also af- 
fected by sums of the form shown in Eq. (22 1, we con- 



clude that, for DTQRWs, the proof of the limitations of 
Lq cannot be applied to L 2 , just as it cannot be applied 
to L\. Therefore, for DTQRWs with L and L\ graph 
certificates, it is possible that there exists a fixed p such 
that the p-particlc nonintcracting walk can distinguish 
all SRGs. 

We will now show that, using the same method of proof 
from the previous subsection, no CTQRW of a fixed num- 
ber of noninteracting particles can distinguish all SRGs, 
when used in conjunction either the L± or L 2 graph cer- 
tificates. We will only examine bosonic walks; the results 
for fermionic walks can be proved in the same manner. 

We begin by stating the definition of L\ for p-boson 
CTQRWs. There are multiple ways in which we can ex- 
tend L\ for more than two particles in a CTQRW, as 
we can choose the number of particles in the final state 
whose probability distribution we wish to measure. For 
simplicity, here we sum over one particle. However, our 
proof holds if we sum over more than one particle as well. 

We examine the certificate L 1 ' p (t): 



where a, (3, and 7 are all functions of the family parame- 
ters and the time t. (The fermionic walks have the same 
form, except the states are anti-symmetrized, and a, j3, 
and 7 are replaced by their respective complex conju- 
gates.) 

Therefore, b (ii ■ ■ ■ ip\ Uc, P (t) \jx ■ ■ -jp) B is a function 
of family parameters, t, and the binary relationships A xy 
and S xy for x G {k . . . i p } and y € {ji ■ ■ - jp}- Moreover, 
this quantity is a sum of terms where each term is a prod- 
uct of up to p binary relationships and various powers of 
a, /3, 7. Each term does not contain any more than one 
instance of any one of the 2p indices. 

Hence \b {k • ■ ■ i p \ \Jc,p(t) \ji ■ ■ ■ jp) B | 2 contains no 
term which includes more than two instances of the 
same index. Thus, J^n \b (k ■ ■ ■ i P \ Uc,p(*) |ji ■ ■■j P ) B I 2 
is only a function of family parameters and the bi- 
nary relationships A xy and S xy for x £ {i 2 ■ ■ ■ i P } 
and y € {ji---jp}- For this quantity to not be 
a function of family parameters, the index k would 
have to appear three times in a single term in 
\ B (k ■■■i P \ U c , P (t) |ji • ■ ■ j P ) B | 2 , as explained in [15]. 

The value of each unique element in Lf' p {t) is de- 
termined by relationships between the 2p — 1 indices 
{12 ■ ■ ■ i p , ji ■ ■ - jp}- Following |15j . it can be shown that 
the number of unique elements L^' p (t) can contain is 
bounded above by 2P 2 +°(p 1o sp). Additionally, the length 



JY 



L c ; > p (t) = sort{]T \b (k... i P \ Uc, P (t) \h 



■h) b 



oiL^ p {t) is { +p )( p p i)- Following the same argu- 
ment from the end of the previous subsection and from 
Section IV C of [TS], the number of unique graph certifi- 
cates this process can generate for a fixed particle number 
p is polynomial in graph size N. Thus a CTQRW used 
with the L\ method and a fixed number of noninteract- 
ing particles cannot be universal for SRGs, as there will 
exist more SRGs than unique L x ' p (t) certificates. 

This analysis extends to L 2 ' P (T). The possible values 
elements of L^' P (T) will in general be different from the 
possible values of Lf' p (t). However, the number of pos- 
sible values elements of L 2 ' P (T) can take on will still 
be bounded above by 2P 2+0 ^ lo ^\ Additionally, the 
lengths of the L^ P {T) and L^ p (t) certificates are the 
same. Thus the number of SRGs of size N in a single 
family that L^' P {T) certificates can distinguish will be 
polynomial in N. We conclude that computing L 2 ' P (T) 
certificates for a fixed number p of noninteracting parti- 
cles cannot distinguish all SRGs. 



1 < ii < is ■ ■ ■ < i P < N, 1 < j x < 31 ■ ■ ■ < j P < N}. 

(23) 

It has been shown [TU [T5] that for SRGs, elements of 
the noninteracting p-boson CTQRW evolution operator 
have the following simple form: 

B (ii...ip|U c , p (£) \j 1 ...j p ) B = 
B (h . . . i v I {at + f3 J + 7 A)®" I jx . . . Jp ) B (24) 



VII. SUMMARY 

We have shown how single-particle discrete-time quan- 
tum random walks can distinguish many non-isomorphic 
strongly regular graphs. These results are proven with 
techniques used to analyze the distinguishing power of 
nonintcracting continuous-time quantum random walks 
[TBTTTS] . These results regarding single-particle DTQRWs 
in turn explain the results of Berry and Wang |18j , who 
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numerically found that two-particle discrete-time walks 
could distinguish many strongly regular graphs. 

Additionally, we have examined a proposal of |18j . 
that DTQRWs have more distinguishing power than 
CTQRWs. To evaluate this proposal, we have found that 
it is important to consider not just the kind of QRW in 
question, but the method in which the graph certificate 
is constructed. We have considered three related graph 
certificate construction methods, which we have dubbed 
L Q , Li, and L 2 . 

We have found that single-particle DTQRWs used with 
the Lq method can distinguish many SRGs, in contrast to 
single-particle and noninteracting two-particle CTQRWs, 
which, when used with the Lq method, cannot distinguish 
any SRGs from the same family, as proven in |13l 114] . 
However, we have also extended the results of 15J, which 
showed that there does not exist a fixed particle num- 
ber p such that a noninteracting p-particle CTQRW with 
the Lo can distinguish all SRGs. Here we have shown 
this limitation to hold true for noninteracting p-particle 
DTQRWs as well. 

Lastly, we have shown that this limitation holds for 
CTQRWs when the L\ and L 2 certificate methods are 
considered. There does not exist a fixed number p such 
that a noninteracting p-particle CTQRW with either L\ 
or L 2 certificates can distinguish all SRGs. However, it 
remains an open question as to whether or not these limi- 
tations of L\ and L 2 apply to DTQRWs. Thus it is possi- 
ble that there exists a noninteracting p-particle DTQRW 
such that Li or L 2 certificates can distinguish all SRGs. 
This would demonstrate a nontrivial difference in dis- 
tinguishing power between continuous-time and discrete- 
time noninteracting walks. 
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IX. APPENDIX 



Eq. (22| in Section 



IV B 



- demonstrates that L®' 1 certifi- 
cates can distinguish non-isomorphic SRGs of the same 
family. To show that Eq. ( 22 1 is of the correct form and 



cannot be a function of family parameters, we provide 



here the value of Ylf=i I Uz?,i(3) \kl) 



(ij\V D ,i(3)\kl) 



N 
N 

^AijAudAik 



(25) 



Su (Jr 



2 s jk 



16 x. 



64/12 ^ 



(^>(A lk + A 3l ) 

I \ A jk) ) 

We recall the SRG identity of Eq. ^ , and see that 
(d — fi)Sjk + p + (A — fj^Ajk- Thus by inspection 



Ojk 



A 2 
A jk 



of Eq. (25), we see that - =1 I U(3) \kl) | 2 contains 



a term proportional to Am Ylf=i AijAjiAjk, corroborat 



ing Eq. (22 1. Therefore, no element of if' 1 can be a 



function only of SRG family parameters, as explained in 
Section MM 
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